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ABSTRACT 

^  Statistics  on  the  frequency  of  transition  between  active  region  classes  were  calculated  and  used  to  derive  the  expected  rate 
of  flaring  for  the  next  2-7  days.  Five  years  of  McIntosh  active  region  classifications  were  analyzed  using  the 
newly-developed  software  called  TELSAR  (Tracking  and  EvoLution  of  Solar  Active  Regions).  The  most  frequent 
transitions  between  these  63  different  classes  usually  involve  a  single-step  change  in  one  of  the  three  classification 
parameters.  The  evolution  of  the  classes  and  the  average,  single-day  flare  rate  were  used  to  predict  flare  rates  for  each 
class  over  the  next  several  days.  For  flares  of  X-ray  class  M  (peak  1~6  X  flux-between  ifl*^^dTO~^~V'Tm~*^these  rates 
can  differ  from  persistence  predictions  by  more  than  0.5  flares  per  day,  although  a  more  typical  difference  is  0. 1  flares  per 
day.  The  rarity  of  some  classes  suggests  that  merging  these  classes  with  others  may  improve  flare  statistics  and  thereby 
improve  flare  forecasts,  f  _ 

1.  INTRODUCTION 

Operational  solar  flare  predictions  are  routinely  made  using  the  reported  three-parameter  McIntosh  active  region  classes 
and  the  flare  probability  for  each  of  these  classes.  A  total  solar  flare  forecast  for  each  day  is  derived  from  the  sum  of  the 
average  number  of  flares  for  the  class  of  each  observed  solar  active  region.  Flare  predictions  for  more  than  one  day, 
however,  require  estimates  of  the  evolution  of  the  region  and  how  much  the  flaring  rate  should  change.  At  present  these 
estimates  rely  on  the  forecaster’s  experience,  because  no  evolutionary  statistics  were  previously  available.  The  present 
study  provides  these  evolutionary  statistics,  and  combines  them  with  flare  probabilities  to  give  an  expected  flare 
probability  for  the  next  2-7  days. 

.Although  every  solar  active  region  exhibits  different  white  light  spatial  characteristics,  each  can  be  assigned  to  a  class, 
using  the  three-parameter  extension  (McIntosh  1981, 1989)  to  the  original  single-parameter  Zurich  sunspot  classification 
(Kiepenheuer  1953,  Waldmeier  1955).  These  three  parameters,  summarized  in  Figure  1,  are  based  primarily  on  the 
region’s  magnetic  configuration,  degree  of  penumbra!  development,  spatial  scales,  and  spot  distribution.  The  first 
parameter  describes  the  nature  of  the  group,  the  second  parameter  describes  the  degree  of  development  of  the  largest 
spot  in  the  group,  and  the  third  parameter  describes  the  amount  of  spot  development  in  the  interior  of  the  group. 
Together,  these  three  parameters  form  a  set  of  63  allowable  classes  (Kildahl  1980),  that  can  be  used  to  examine  the 
complexity  and  maturity  of  the  region.  We  have  listed  these  parameters  from  left  to  right  and  then  top  to  bottom  in  the 
summary  figure  in  rough  order  of  increasing  complexity.  This  is  not  a  definitive  ordering,  and  future  analysis  may  suggest 
reordering  these  parameters. 

Because  of  solar  rotation,  the  entire  evolution  of  a  single  solar  active  region  is  rarely  observable.  We  assume  that 
similarities  of  evolution  occur,  and  that  the  evolution  of  each  region  is  not  totally  random.  This  assumption  is  based  on  the 
evolutionary  nature  of  the  original  single-parameter  classification  (Bray  and  Loughhead  1964).  Therefore,  by  comparing 
the  evolution  of  many  regions,  we  can  derive  the  probable  evolution  of  a  solar  active  region. 

2.  THE  TELSAR  ANALYSIS  SOFTWARE 

We  have  recently  developed  new  software,  called  TELSAR  for  Tracking  and  EvoLution  of  Solar  Active  Regions,  to 
calculate  statistics  on  the  evolution  of  active  regions  through  the  three-parameter  classification  system.  For  our  analysis. 
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Figure  1.  Schematic  representation  of 
criteria  used  to  define  each  of  the  three 
McIntosh  active  region  parameters. 
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we  have  added  two  more  classifications  to  examine  the  evolution  from  birth  to  death  of  a  region.  The  birth  class,  YYY, 
was  added  if  the  region  was  first  observed  west  of  E70  deg.  Similarly,  a  death  class,  ZZZ,  was  added  if  the  last  region  report 
occurred  before  reaching  W70  deg. 

The  reported  active  region  classes  for  each  region  were  converted  into  a  separate  evolutionary  “tree”  showing  the 
evolution  from  one  class  to  the  next.  These  trees  were  then  used  to  calculate  evolutionary  and  flare  statistics.  The 
evolutionary  statistics  include  the  minimum,  maximum,  and  average  length  of  time  for  regions  of  one  class  to  evolve  into 
another  class,  and  probability  and  total  number  of  regions  making  these  transitions.  The  flare  statistics  give  the  probable 
number  of  flares  per  day  for  each  class  after  a  specified  number  of  days  have  elapsed. 

Because  a  region  was  usually  classified  by  a  number  of  different  observatories,  the  classifications  for  a  given  day  may  be 
inconsistent.  TELSAR  has  two  primary  methods  for  dealing  with  these  inconsistencies.  The  first  method  is  an  >0(j#a 
equal-weight  vote  from  each  of  the  reported  classifications  on  that  day.  A  tie  vote  is  resolved  using  the  reported  reliability  ^  Qr 


of  the  observations.  The  second  method  is  a  reliability-weighted  average  of  each  of  the  three  classification  parameters. 
Combinations  of  parameters  that  would  give  an  illegal  class  are  redefined  according  to  an  algorithm  used  by  the  SEL 
Space  Environment  Services  Center  (developed  by  J.  Bryson  and  D.  Snelling;  C.  Cruickshank  1988  private 
communication).  We  have  used  this  second  method,  the  SESC  average,  for  all  of  the  analysis  presented  in  this  paper. 

I'ELSAR  was  designed  to  analyze  subsets  of  the  data  based  on  different  criteria.  Data  can  be  chosen  based  on  a  range  of 
dates,  range  of  region  numbers,  or  evolution  through  a  particular  class.  Less  reliable  data  can  be  automatically  discarded 
whenever  the  observatory-reported  reliability  lies  below  some  limit.  Data  can  also  be  discarded  whenever  the  region  was 
within  a  specified  distance  from  the  limb,  where  foreshortening  makes  the  classification  more  difficult. 

TELSAR  can  compare  the  evolutionary'  patterns  of  the  separate  trees  and  combine  them  when  they  have  a  sufficient 
number  of  similar  transitions  between  classes.  This  is  done  tn  an  effort  to  piece  together  the  probable  evolution  of  regions 
on  the  back  side  of  the  sun.  The  resulting  combined  trees  (shown  schematically  in  Figure  2)  can  then  be  analyzed  to 
determine  the  statistical  properties  beyond  the  observational  window. 
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Because  TELSAR’s  ability  to  correctly  combine  evolutionary  sequences  is  still  being  tested,  this  paper  will  concentrate  on 
statistics  derived  without  combining  trees.  Because  no  tree  combining  is  performed,  the  evolution  into  the  far  future  will 
not  be  considered.  The  evolution  several  days  into  the  future  should  not  be  severely  affected,  since  the  observational 
window  due  to  solar  rotation  is  approximately  14  days. 

3.  DATA 

The  present  analysis  uses  a  subset  of  the  active  region  classifications  available  on  magnetic  tape  from  the  World  Data 
Center.  We  have  used  the  classifications  from  December  1981  through  the  end  of  1986,  which  includes  a  total  of  1158 
active  regions  with  identification  numbers  between  3477  and  4762.  A  total  of  32,932  individual  active  region  reports  were 
examined,  of  which  32,705  contained  useable,  legal  classifications  (see  section  4.3).  The  three-parameter  active  region 
classifications  were  assigned  in  real  time  at  the  U.S.  Air  Force  Solar  Observing  Optical  Network  (SOON)  sites  for  use  in 
flare  forecasting  by  the  NOAA  Space  Environment  Services  Center  and  the  USAF  Global  Weather  Central.  Up  to  eight 
different  sites  contributed  their  daily  active  region  reports,  which,  after  averaging  all  reports  for  the  same  region  and  day, 
gave  a  total  of  7030  region  classifications.  Therefore,  an  average  of  4.7  daily  reports  were  combined  to  form  an  average 
sequence  of  6. 1  days  of  classifications.  We  have  chosen  not  to  analyze  the  early  10  years  of  data  at  this  time,  in  part  because 
many  of  the  SOON  sites  were  not  yet  operational. 

4.  RESULTS 

4.1  TRANSITIONS  BETWEEN  CLASSES 

In  Table  1  we  present  the  statistics  on  the  number  of  direct  transitions  between  classes.  This  table  is  a  65x65  array  listing 
the  scaled  number  of  regions  that  evolved  from  the  class  listed  in  the  left  hand  side  into  the  class  listed  at  the  top.  We  have 
ordered  the  classes  from  the  left  to  right  and  then  top  to  bottom  as  shown  in  Figure  1.  Appended  to  the  end  of  this  class 
sequence  are  our  birth  and  death  classes.  Each  class  was  assigned  a  class  number,  given  in  parenthesis.  These  same 
two-digit  class  numbers  are  repeated  at  the  top  of  the  table,  and  should  be  top  to  bottom. 

In  order  to  conserve  space,  we  have  scaled  the  number  of  direct  transitions  between  SESC-avcragcd  classes  and  rounded 
this  scaling  to  a  integer.  We  use  a  logarithmic  scaling,  s,  of  the  form 

s  =  10.  *  log  x  /  log  [max(x)| 

where  x  is  the  number  of  regions  observed  to  make  a  direct  transition  between  classes.  The  scaled  numbers  run  from  0  to  9 
with  a  number  sign  for  the  maximum  of  10.  The  maximum  number  of  transitions,  max(x),  was  431,  so  this  scaling 
reduces  to 
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X  =  1.83s 

and  the  scaled  numbers  from  0  to  4  correspond  to  approximately  1,  2, 3,  6,  and  1 1  occurrences.  In  addition,  a  period  is 
used  to  indicate  transitions  that  did  not  occur,  and  an  asterisk  notes  the  classes  that  did  not  occur  in  the  original  data. 

A  large  number  of  transitions  occurred  between  neighboring  classes.  For  ~  80%  of  the  classes,  the  most  likely  next  class  is 
only  one  parameter  different.  And  in  these  cases,  the  one  parameter  generally  changed  by  only  one  step  in  the 
left-to-nght  progression  of  the  flare  parameters  shown  in  Figure  1.  The  primary  exception  to  this  is  the  set  of  second 
parameters  S.  H,  A,  and  K,  which  are  based  on  combinations  of  the  sunspot  size  and  symmetry. 

For  one  third  of  the  classes,  the  most  probable  transition  into  and  out  of  this  class  involved  the  same  second  class.  This 
suggests  the  two  possibilities  illustrated  in  Figure  3.  Either  the  reported  classes  are  oscillating  between  neighboring 
classifications  due  to  observational  errors,  or  the  sequence  followed  by  growing  regions  is  traversed  in  the  opposite 
direction  by  decaying  regions.  These  two  possibilities  can  be  tested  using  a  running  average  of  the  classes  with  time.  Such 
an  average  should  reduce  the  number  of  oscillating  classifications.  Alternatively,  the  examination  of  individual  or 
combined  trees  may  reveal  the  answer.  These  tests  have  not  yet  been  performed. 


Class  1  Class  2 


Figure  3.  Traversing  the  same  transitions  forwards  and  then  backwards  or  oscillating  between 
classes  are  two  possible  causes  for  the  result  that  the  most  probable  transitions  into  and  out 
of  a  class  often  involve  the  same  intermediate  class. 

4.2  FLARE  PROBABILITIES 

I'ELSAR  calculates  the  rate  of  M  flares  per  day  expected  to  occur  2  to  7  days  after  the  observed  classification.  These  rates 
were  derived  using  all  observed  active  regions  from  December  1981  through  1986,  without  the  combining  indicated  in 
Figure  2.  For  each  class,  TELSAR  determined  all  classes  that  those  regions  evolved  into  after  a  specified  number  of  days. 
Ibis  was  then  used  along  with  the  average  number  of  flares,  P  j,  for  these  subsequent  classes  (Kildahl  1980)  to  calculate 
the  expected  flare  rate  for  that  day.  (We  used  his  table  of  class  occurrences  to  resolve  two  typographical  errors  in  the  flare 
probability,  which  listed  the  flare  probabilities  for  classes  HRX  and  FRO  listed  twice,  and  omitted  HHX  and  FKO.)  Thus, 
this  is  not  the  actual  number  of  flares  observed  for  that  region  several  days  later,  but  a  statistically  predicted  average. 
Finally,  the  expected  flare  rate,  R,  in  units  of  M  flares  per  day  is  derived  from 

R  =  Z  Pj  Nj  /£  Nj 

were  N  j  is  the  number  of  regions  of  class  i  observed  on  the  specified  day.  We  have  included  a  zero  average  number  of 
l  lares  for  the  birth  and  death  classes,  but  did  not  consider  a  region  to  contribute  zero  flares  after  t*e  first  day  in  the  death 
class  (ZZZ).  Ihis  could  affect  the  calculated  flare  rate,  and  cause  an  artificial  rise  for  the  day  after  the  death  class. 
1  lowever,  smee  decaying  regions  do  not  produce  many  flares,  this  will  not  seriously  alter  the  Hare  forecast.  We  also  note 
that  the  confidence  in  this  flare  rate  will  decrease  with  time  because  the  number  of  regie  .s  observed  for  long  periods  of 
time  decreases  due  to  solar  rotation  with  the  time  interval. 

In  table  2  we  list  the  calculated  flare  rates.  The  first  column  lists  the  classes,  along  with  an  asterisk  for  those  classes  that  did 
not  occur  in  the  original  data.  The  second  column  lists  the  average  number  of  flares  per  day  derived  by  Kildahl  (1980)  for 
the  years  1969  to  1976.  The  remaining  columns  list  the  future  flare  rate  R,  wit  ii  periods  to  indicate  when  no  classes  were 
observed  on  the  specified  future  day. 

Regions  near  the  top  of  the  table  tend  to  be  simple,  and  produce  increasing  numbers  of  flares  with  time.  Regions  near  the 
bottom  of  the  table  arc  more  complex  and  are  likely  to  have  a  maximum  flare  rate  within  the  next  several  days. 

These  evolutionary  statistics  are  expected  to  provide  better  flare  forecasts  for  the  next  several  days  than  persistence 
forecasts.  The  presistence  forecast  does  not  consider  evolution  of  the  region  or  changes  in  the  flare  rate  with  time.  This 
can  be  seen  by  comparing  the  same  day  flare  rates  in  column  two  ofThble  2  with  the  following  columns.  The  maximum 
flare  rate  expected  in  the  next  several  days  can  be  up  to  a  factor  of  7  larger  than  for  the  current  day.  although  a  more  typeial 
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Table  2.  Average  number  of  M  flares  per  day 


Class 

Per  dav 

day  ,2 

day  3 

dav  4 

day  5 

day  6 

AXX (  0 ) 

0.01 

0.02 

0.03 

0.04 

0.06 

0.08 

HRX (  1) 

0.03 

0.03 

0.04 

0.03 

0.06 

0.07 

HSX (  2) 

0.05 

0.06 

0.06 

0.07 

0.08 

0.09 

HHX (  3) 

0.11 

0.09 

0.12 

0.13 

0.12 

0. 13 

HAX (  4) 

0.06 

0.06 

0.11 

0.11 

0.12 

0. 11 

HKX (  5) 

0.18 

0.20 

0.22 

0.24 

0.20 

0.22 

BXO (  6) 

0.02 

0.03 

0.03 

0.04 

0.06 

0.08 

BXI (  7) 

0.06 

0.12 

0.40 

0.26 

0.45 

0.15 

CRO (  8) 

0.05 

0.04 

0.05 

0.06 

0.08 

0.08 

CRI (  9) 

0.05 

0.17 

0.27 

0.18 

0.23 

0. 32 

CSO(IO) 

0.04 

0.06 

0.06 

0.06 

0.07 

0.07 

CSI(ll) 

0.08 

0.10 

0.11 

0.06 

0.03 

0.02 

CHO ( 12 ) 

0.07 

0.09 

0.12 

0.14 

0.13 

0.13 

CHI (13) 

0.21 

0.18 

0.15 

0.09 

0.19 

CAO ( 14 ) 

0.08 

0.08 

0.09 

0.12 

0.12 

0.12 

CAI(15) 

0.11 

0.08 

0.10 

0.04 

0.01 

CKO (16) 

0.25 

0.16 

0.19 

0.15 

0.14 

0.15 

CKI (17) 

0.21 

0.31 

0.28 

0.37 

0.10 

0 . 49 

DRO ( 18 ) 

0.08 

0.05 

0.04 

0.07 

0.04 

0.03 

DRI (19) 

0.13 

0.10 

0.18 

0.08 

0.05 

0.00 

DRC ( 20) * 

0.00 

DS0(21) 

0.09 

0.10 

0.12 

0.12 

0.14 

0. 16 

DSI (22) 

0.13 

0.21 

0.25 

0.24 

0.22 

0.29 

DSC (23) 

0.25 

DHO ( 24 ) 

0.26 

0.16 

0.27 

0.24 

0.21 

0.20 

DHI ( 25 ) 

0.07 

0.48 

0.37 

0.36 

0.41 

0.35 

DHC ( 26 ) 

0.33 

DAO (27) 

0.10 

0.14 

6.16 

0.17 

0.17 

0.16 

DAI (28) 

0.18 

0.30 

0.31 

0.32 

0.33 

0.29 

DAC(29) 

0.26 

0.08 

0.08 

0.10 

0.10 

0.08 

DK0(30) 

0.33 

0.31 

0.30 

0.34 

0.36 

0.24 

DKI (31) 

0.48 

0.49 

0.53 

0.46 

0.49 

0.45 

DKC ( 32 ) 

0.72 

0.66 

0.55 

0.50 

0.50 

0.45 

ERO (33) 

0.00 

0.09 

0.05 

0.11 

0.10 

0.08 

ERI (34) 

0.00 

0.08 

0.05 

0.05 

0.02 

0.00 

ERC ( 35 ) * 

0.00 

ESO(36) 

0.17 

0.16 

6',  i3 

0.17 

0.13 

0.18 

ESI (37) 

0.17 

0.24 

0.26 

0.31 

0.27 

0.16 

ESC (38) * 

0.00 

EHO (39) 

0.15 

0.26 

0.38 

0.34 

0.34 

0.29 

EHI (40) 

0.62 

0.33 

0.35 

0.23 

0.16 

0.19 

EHC (41) 

2.00 

EAO (42) 

0.21 

0.27 

0.25 

0.30 

0.33 

0.23 

EAI (43) 

0.58 

0.58 

0.53 

0.51 

0.47 

0. 48 

EAC(44) 

0.35 

0.35 

1.27 

0.58 

0.10 

0 . 00 

EKO(45) 

0.38 

0.40 

0.50 

0.43 

0.36 

0 . 44 

EKI (46) 

1.27 

1.02 

0.96 

0.89 

0.81 

0 . 78 

EKC (47) 

2.37 

2.07 

1.71 

1 . 50 

1.16 

1.09 

ERO (48) 

0.00 

FRI (49) * 

0.50 

FRC(50) * 

0.00 

FS0(51 ) 

0.46 

0.23 

0.40 

6.27 

0.10 

0.36 

FSI(52) 

1.88 

0.00 

0.00 

0.09 

0.06 

FSC (53 ) * 

0.00 

FHO(54) 

0.00 

0. 1.4 

0.09 

6.21 

0.29 

0. 16 

EHI (55) 

0.83 

1.01 

0.50 

0.65 

0.89 

0.30 

EHC (56 ) 

0.80 

E  AO (57) 

0.00 

0.1.5 

0.34 

0.07 

0.06 

0.08 

EAI (58) 

0.67 

1.23 

0.55 

0.62 

0.15 

0.05 

FAC (59) 

0.00 

EKO ( 60 ) 

0.32 

0.79 

0.65 

0.44 

0.55 

0. 62 

FKI (61 ) 

2.26 

1.64 

1 . 49 

1.21 

1.00 

0.98 

EKC(62) 

1.44 

1.85 

1.92 

1.37 

1.19 

1 . 01 

YYY ( 63 ) 

0.00 

0.03 

0.04 

0.05 

0.06 

0.08 

ZZZ ( 64 ) 

0.00 

day  7 
0.07 
0.14 
0.11 
0.15 
0.08 
0.26 
0.09 
0.30 
0.10 
0.12 
0.09 
0.02 
0.12 

0.11 

0.12 

0.34 

0.04 


0. 17 
0.31 

0.19 

0.18 

0.17 

0.29 

0.08 

0.32 

0.46 

0.45 

0.08 


0.22 

0.08 

0.42 

0.15 

0.28 

0.43 

0.51 

0.66 

1.16 


0.10 


0.23 

0.55 

0.08 

0.20 

0.75 

0.87 

0.82 

0.10 
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differences  is  a  factor  of  2.  Alternatively,  the  maximum  flare  rate  over  the  next  several  days  typically  differs  from  the 
persistence  rate  by  0.1  flares  per  day,  or  in  more  extreme  cases  by  0.6  flares  per  days. 

We  note,  however,  that  the  most  extreme  differences  between  these  statistics  and  persistence  predictions  generally  occur 
for  rarely  observed  classes.  These  differences  may  be  an  artifact  of  poor  statistics  for  these  classes.  This  would  suggest 
that  merging  these  rare  classes  with  similar  classes  may  improve  the  ability  to  acurately  predict  flares  from  these  regions. 

4.3  ILLEGAL  AND  RARE  CLASSES 

Analysis  of  the  data  pointed  to  several  types  of  infrequent  errors  in  the  reported  active  region  classes.  The  most  common 
tvpe  of  error  were  illegal  classes.  In  these  cases,  summarized  in  Table  3.  two  of  the  parameters  disagreed  about  either  the 
existence  of  penumbra  or  the  polarity  of  the  region.  These  illegal  classes  occurred  201  times,  which  is  only  0.6%  of  all  of 


the  reported  classes. 

Table  3.  Incompatible 

Class  Parameters 

Parameter 

Unipolar  Bipolar 

No  Penumbra 

Penumbra 

first 

AH  B.  C.  D,  E.  F 

A  B 

H,  C,  D,  E,  F 

second 

X 

R,  S,  H.  A  K 

third 

X  O.  I,  C 

I.  C 

The  second  type  of  classification  errors  were  two  classes.  CKC  and  CSC,  that  were  not  clearly  illegal  based  on  the  early 
definitions  of  these  classes  (McIntosh  1986  and  the  guidelines  for  observers  in  the  Air  Weather  Servce  Pamphlet  105-60). 
The  classes  are  now  listed  as  illegal  in  the  classification  paper  of  McIntosh  (1989),  probably  because  it  is  somewhat 
contradictory  to  require  penumbra  on  spots  in  the  middle  and  at  only  one  end  of  the  region.  We  have  treated  the  four 
occurrences  of  these  two  classes  as  illegal,  and  did  not  process  these  reports. 


A  third  type  of  error  were  classes  in  which  one  or  more  of  the  three  classification  parameters  were  missing.  Of  these, 
there  were  22  cases  that  could  not  be  completed  from  the  existing  parameters.  The  incomplete  classes  could  be  resolved 
from  the  two  existing  parameters,  using  the  criteria  in  "Ihble  3.  These  classes  were  corrected  to  the  full  three-parameter 
class  (see  Table  4),  and  the  analysis  of  these  regions  continued  as  normal. 


Table  4.  Resolution  of  Incomplete  Classes 


B.saQr..t£d 

AX_ 

B_0 

b_i 

hr_ 

HS_ 

HH_ 

HA_ 

HK 


Complete 

AXX 

BXO 

BXI 

HRX 

HSX 

HHX 

HAX 

HKX 


A  fourth  type  of  error,  which  is  less  easily  recognized  is  the  improper  assignment  of  a  region  to  a  legal  class.  These 
improper  classifications  can  occur  due  to  poor  seeing  or  observing  conditions,  foreshortening  of  regions  on  the  limb,  and 
human  error  (for  further  discussion  see  McIntosh  1989).  Because  the  forecaster  must  make  predictions  based  on  data  that 
may  contain  these  errors,  we  felt  it  best  to  calculate  our  statistics  on  this  “noisy"  data.  Therefore  we  have  not  attempted  to 
identify  these  errors. 

Twelve  classes  did  not  occur  in  the  SESC-average  of  the  1981-1986  data  (see  rows  and  columns  of  periods  in  Table  1).  We 
searched  the  entire  original  data  set  for  the  number  of  occurrences  of  these  twelve  classes.  This  extended  data  set  included 
all  of  1971  through  1986,  and  was  not  limited  to  the  smaller  sample  considered  in  the  rest  of  this  paper.  These  classes  are 
listed  in  Table  5  in  a  way  to  suggest  combinations  of  parameters  that  should  possibly  be  considered  illegal.  One  such 
combination  has  FR  for  the  first  two  parameters,  which  requires  a  large  penumbra  for  the  major  spot  and  only 
rudimentary  penumbra  for  the  interior  spots.  Another  set  of  combinations  have  RC,  SC,  or  HC  as  the  last  two  parameters, 
which  requires  cither  rudimentary  or  small,  symmetric,  mature  penumbra  for  the  major  spot  and  mature  penumbra  for 
the  interior  spots.  We  note  that  the  RC  combination  is  now  considered  illegal  by  McIntosh  (1989),  although  nodetailed 
reasons  were  presented.  We  suggest  that  these  different  combinations,  along  with  possibly  class  FAC.  may  indicate  a 
conflict  about  the  degree  of  complexity  and  development  of  the  region. 
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FAC  7 


Thble  5.  Rare  classes  and  number  of  occurrences  in  1971-1986 

FRO  4  DRC  1  DSC  14  DHC  12 

FRI  2  ERC  0  ESC  1  EHC  5 

FRC  0  FSC  0  FHC  3 

Because  of  the  rarity  of  these  classes,  a  single  improperly  assigned  classification  or  improperly  assigned  region  responsible 
for  the  flare  could  increase  the  average  number  of  M  flares  by  up  to  0.5  flares  per  day.  Merging  these  rare  classes  with 
similar  classes  would  reduce  the  formal  errors.  If  the  classes  are  sufficiently  similar  that  they  could  easily  be  confused,  or 
might  rapidly  evolved  from  one  to  the  other,  then  such  a  merger  would  not  introduce  new  errors  into  the  flare  rate  of  the 
more  frequent  class. 

5.  CONCLUSIONS 

We  have  presented  the  initial  results  of  the  newly-developed  analysis  software,  TELSAR.  This  software  has  been  applied 
to  5  years  of  active  region  classifications  to  derive  statistics  on  the  frequency  of  transition  between  classes  and  the 
subsequent  rate  of  flaring  expected  from  the  region  for  the  next  2-7  days.  These  tabulated  statistics  form  the  first  study  of 
this  kind  that  we  are  aware  jf.  It  is  hoped  that  these  results  can  be  used  to  improve  flare  forecasts  on  the  intermediate 
timescale  of  several  days.  ' 

The  most  frequent  transitions  between  classes  involve  single-step  changes  in  only  one  of  the  three  classification 
parameters.  The  reverse  of  these  transitions  were  also  very  frequent,  suggesting  that  the  regions  either  oscillate  in  class 
due  to  observational  error,  or  traversed  the  same  sequence  forwards  and  then  backwards  as  the  region  grows  and  then 
decays.  This  could  be  tested  with  temporally  averaged  classes.  In  addition,  such  averages  may  reduce  the  number  of 
entries  in  the  transition  matrix  by  reducing  the  frequency  of  rarely  observed,  spurious  transitions. 

The  maximum  future  flare  rate  differs  from  persistence  predictions  by  up  to  0.5  M  flares  per  day,  with  0.1  being  a  more 
typical  difference.  These  initial  results  suggest  that  evolutionary  statistics  can  be  used  to  improve  intermediate  term  flare 
forecasts,  but  their  significance  has  not  yet  been  tested.  One  such  test  would  compare  the  number  of  flares  actually 
produced  by  the  region  several  days  later  with  the  TELSAR-derived,  expected  rate.  Another  test  would  determine  the 
difference  between  the  TELSAR-derived  future  rates  and  the  persistence  rates.  This  could  be  done  through  an  analysis  of 
the  propagation  of  errors,  starting  with  the  statistical  errors  on  the  average  number  of  flares  observed  per  day  for  each 
class,  i  he  resulting  errors  on  the  TELSAR-derived  flare  rates  should  exceed  the  difference  with  persistence  rales  for  the 
evolutionary  statistics  to  prove  valuable. 

The  frequency  of  observation  of  several  classes  was  very  low,  often  less  frequent  than  reports  of  illegal  classes.  This 
suggests  that  reduction  in  the  number  of  allowed  classes  should  be  considered.  Such  a  reduction  may  improve  the  flare 
statistics  and  thereby  improve  flare  forecasts. 
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